Analytical approach to soliton ratchets in asymmetric potentials 
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We use soliton perturbation theory and collective coordinate ansatz to investigate the mechanism 
of soliton ratchets in a driven and damped asymmetric double sine-Gordon equation. We show that, 
at the second order of the perturbation scheme, the soliton internal vibrations can couple effectively, 
in presence of damping, to the motion of the center of mass, giving rise to transport. An analytical 
expression for the mean velocity of the soliton is derived. The results of our analysis confirm the 
internal mode mechanism of soliton ratchets proposed in [Phys. Rev. E 65 025602 (R) (2002)]. 
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I. INTRODUCTION 

During the past years a great deal of attention has 
been devoted to the ratchet effects both for point parti- 
cles nil, nil and for extended systems [IIEH- Some 
experimental realization of these models can be found 
in 0, 1^ iEj ilAl JJJ fo'' ^ recent review). Ratchet- 
like systems such as systems of two particle with internal 
degree of freedom B, periodic rocket ratchets (determin- 
istic and stochastic) ^3| and temperature ratchets [isj, 
have also been considered. For infinite dimensional sys- 
tems described by nonlinear partial differential equations 
(PDE) of soliton type, ratchet phenomena were investi- 
gated both in the case of asy mmetric potentials in pres- 
ence of symmetric forces , and in the case of symmet- 
ric potentials with asymmetric forces 16] . In both 
cases the ratchet phenomenon manifests as an unidirec- 
tional motion of the soliton, similar to the drift motion 
occurring for point particle ratchets (from here the name 
of soliton ratchets). A symmetry approach to the phe- 
nomenon, which allows to establish conditions for the oc- 
currence of soliton ratchets, was developed in Ref. [lij . 
This approach, although useful for predicting the phe- 
nomenon, does not provide information about the ac- 
tual mechanism responsible for the unidirectional mo- 
tion. The mechanism underlying soliton ratchets was re- 
cently proposed by two of us for the case of a perturbed 
asymmetric double sine- Gordon equation driven by sym- 
metric forces Q , and extended in Ref. jl^ to the case of 
a perturbed sine Gordon system in presence of asymmet- 
ric drivers. In both cases the phenomenon was ascribed 
to the existence of an internal mode mechanism which, 
in presence of damping, is able to couple soliton inter- 
nal vibrations to the translational mode, producing in 
this way transport. The phenomenon was described as 
follows: the energy pumped by the ac force into the soli- 
ton internal mode is converted into net dc motion by the 
coupling of the internal vibration with the center of mass 
induced by the damping. This internal mode mechanism 



was confirmed by numerical and analytical investigations 
for the case of the sine Gordon system with asymmetric 
periodic forces 1^ .IS . ,11 \ and by direct numerical inves- 
tigations for the case of soliton ratchets of the asymmet- 
ric double sine-Gordon equation In the last case, 
a detailed numerical investigation (see 0]) showed the 
following facts: 

(i) The presence of damping and the asymmetry of the 
potential are both crucial ingredients for the existence of 
the net motion of the soliton (in presence of the ac force, 
but in absence of damping, the asymmetry of the poten- 
tial does not produce transport), (ii) The effect of the 
damping is to couple the internal vibrations of the kink 
to the motion of the center of mass, (iii) For fixed values 
of the amplitude and frequency of the ac force there is an 
optimal value of the damping for which the transport is 
maximal, (iv) The direction of the motion is fixed by the 
asymmetry of the potential and is independent on initial 
conditions, (v) For fixed values of the damping the av- 
erage velocity of the kink shows a resonant behavior as 
a function of the frequency of the ac force . (vi) At low 
damping and higher forcing strengths, current reversals 
in the kink dynamics can occur. This phenomena was 
ascribed in Ref. |3| to the soliton-phonon interaction 
rather than to the internal mode mechanism. 

These points were found to be valid also for different 
soliton ratchet systems 0, 0| . In Ref. , the exis- 
tence of an optimal value of the damping constant which 
maximizes the transport, was ascribed to the relativis- 
tic nature of the kink dynamics and the possibility of 
nonzero currents in absence of dam ping , was also found. 
Since most of the analysis of Ref. [l^ has been based 
on numerical simulations, it is of interest to provide an 
analytical confirmation of the above results. 

The aim of this paper is to present an analytical in- 
vestigation of soliton ratchets of the asymmetric double 
sine-Gordon equation which confirms the aforementioned 
internal mode mechanism proposed in Ref. (l3| as well 
as points (i) - (vi) listed above. To this end, we use per- 
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turbation theory and a collective coordinate ansatz for 
the soliton shape to derive ordinary differential equations 
(ode) for the center of mass and for the soliton width. We 
show that, to second order of perturbation theory, soli- 
ton width vibrations couple effectively to the motion of 
the center of mass via the damping in the system. By 
using the collective coordinate equations we are able to 
derive an analytical expression of the mean drift velocity 
of the kink as a function of the system parameters. The 
analysis is shown to be in good agreement with numeri- 
cal simulation and with the results (i) - (vi) obtained in 
Ref. 14]. 

The paper is organized as follow: in section II we in- 
troduce our model and discuss its main properties. In 
section III we derive the dynamical equations for soliton 
collective coordinates and obtain an expression of the 
average kink velocity as a function of the system param- 
eters. In section IV we compare our analytical results 
with direct numerical simulations and discuss them in 
connection with previous work. Finally, in section V the 
main conclusions of the paper are summarized. 



II. THE MODEL 

Let us consider the following perturbed asymmetric 
double sine-Gordon equation (ADSGE) 

(l>tt-(l)xx + sin((/)) Acos(20) = F(a;,i, 0, 0t, ...) EE 
fit) - a4>t, (1) 

where A £ [—1,1] is a parameter related to the asymme- 
try of the nonlinear Klein-Gordon potential, a is a damp- 
ing constant and f{t) — esin(ajt -f- ^o) is a periodic force 
with amplitude e, frequency w and phase ^o- This sys- 
tem is connected with interesting physical problems such 
as arrays of inductively coupled asymmetric SQUIDs of 
the type considered in Ref. 19]. A mechanical analogue 
of Eq. ^ in terms of a chain of double pendulum was 
given in Ref. For F = Q, Eq. |^ has an hamiltonian 
structure with Hamiltonian (energy) 
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dx4>x:4>t- (3) 



The potential in Eq. (0) is U{<j)) = C - cos{(j>) + 
(A/2)sin(20) and C = cos((/)o) - (A/2) sin(20o) (notice 
that for A 7^ the potential is asymmetric in the field 
variable). Similar to the sine-Gordon equation, the en- 
ergy and the momentum are both conserved quantities 
for the unperturbed ADSGE. In this case one can show 



[IJ that Eq. 
of the form 



Q has an exact kink (anti-kink) solution 




where ^ — 
77 = 2Av/2(T 
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vt)/VT^ 

B = V'2(4A2" 



1 + A) anu (po 
arcsin[(l — A)/(4A)] -|- 27rn (the plus and minus signs re- 
fer to the kink and antikink solutions, respectively). In 
the limit A ^ (zero asymmetry), Eq. Q reduces to 
the well known soliton solution of sine-Gordon equation 
(SGE). 



III. COLLECTIVE COORDINATE ANALYSIS 

The term F{x,t, 4>, (jjt) added to the ADSGE can be 
considered as a small perturbation acting on the system. 
In this case the energy and the momentum will depend 
on time so that it is natural to assume an ansatz for the 
perturbed kink of the form 



2 tan 



sign(A) ^i? 



A — 1 — rj sinh 



x-X[t) 
W{t) 



(5) 



where X{t) and W{t) represent dynamical collective co- 
ordinates (CC) corresponding to the center of mass and 
the width of the kink, respectively (similar approach was 
introduced in [i^] for the SGE). In the following we 
consider only kink solutions since the analysis of anti- 
kinks will follow from it without difficulty. In absence 
of perturbation the kink moves with constant velocity 
V and co nstant w idth W{t) ^ Ws = Wq^/I - ( 
Wo = 2^\X\/[AB] ), so that X{t) = Vt. Notice that 
ansatz © is consistent with the linear stability analysis 
performed in Ref. |l4| . showing the existence of an in- 
ternal mode frequency 0/ below the phonon band (this 
mode disappears in the sine-Gordon limit A = 0). From 
the spatial profile of the corresponding localized eigen- 
function, one can see that the internal mode is linked to 
vibrations of the kink width, this suggesting the choice 
of ansatz ((SJ. 

By substituting Eq. © into Eqs. l(2Jl, ®, and differ- 
entiating with respect to time, one obtain, after straight- 
forward calculations (for details see H^), that X(t), P(t) 
and W{t) satisfy the following system of nonlinear ode 



dX _ P{t)W{t) 
~dt ~ 
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-aP-qf{t), 



where the momentum P is a solution of the equation 

dP 
dt 

and 



(7) 



[(A - l-7?sinh[6l])2+A2S2]2 



h = 



COS 



[A -I - r;sinh[6i])2 + ^2^2' 



i = l,2,3; 



/i(A)A3g 3^2 



K=\h-'f^\; R = 2ABrj, g ^ 



Notice that the above integrals depend on the asymmetry 
parameter A and are finite for any value of A (the y c an 
be easily calculated by numerical tools). In Ref. [13 a 
similar approach was used to show that the internal mode 
mechanism give rise to net motion in the SG equation 
perturbed by two harmonic forces. 

In order to solve the system of equations (|6a|l -l)6b |l and 
^ we use the following initial conditions: dX{Q)/dt = 0, 
P(0) — and W{0) — Wo- We solve first the non-coupled 
linear equation for the momentum. From Eq. ((TJ we 
obtain 



P{t) = 



-qe 



[F(0) 



y/a^ + Ll) 



sin[wt + 9o + S] + 

sm{9o + 6)l (8) 



where tan((5) — —uj/a and the last term in the r.h.s, being 
a transient, will be neglected (we are interested in the 
stationary regime t >> a~^). Notice from Eq. (|6bl) that 
the width of the kink is affected directly by the ac force 
with a frequency uj and indirectly by the translational 
motion (the momentum) with a frequency 2uj (see Ref. 
[i^b). Since it is difficult to find an exact solution of this 
equation, we will search for an approximate solution for 
W{t) in the form of a power series in e 

W{t) ^Wo + eWi{t) + e^W2{t) + O(e^). (9) 

Inserting into Ij6bp and taking terms of the same order 
in e, we obtain, after a transient time {t » a~^), the 
following set of linear equations. At order 0(e) we have 



Wi + ^fWi + aWi = WqG{X) sm{ujt + 9^), (10) 



where 
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2g(A) 
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qh - sign{X)IihR 



while at the second order (0(e2)) we obtain 



(11) 

(12) 



W2 + n^W2 + aW2 = 



2Wa 2Wa 



(13) 



W„q^ sm^[ujt + 9o + 6] 

2i?4(/i/3 -/2)(a2 +^2- 



+ G{X)Wisiii{ujt + 9o). 



Eqs. (|10|l and H13|l correspond to linear, damped and 
driven oscillators with characteristic frequency fi. It can 
be shown that for A E [— 1 , 1] the values of O are quite 
close to the internal mode frequency fi/ of the ADSGE 
computed in |l4j (the maximum difference between them 
being « 0.1). Eqs. 1(111)1 and l(T^ can be solved exactly. 
The homogeneous part of the solutions of both equations 
decays to zero, after a transient time, if — < 0. 
The particular solutions of lflT)| and (fl^ are given, re- 
spectively, by 



W2{t) 



WoGjX) cos(^t + 9o + 9) 
. ^2 + 

f72[(f]2 _^2)2 +Q,2^2] ^ 

Ax sin(2wi -I- 29^ + 2^ + ^1) 



(14) 
(15) 



[(1^2 _ t^2)2 + Q,2cj2] ^(1^2 _ 4^2)2 +4^2^2 

A2 sin(2wi + 26*0 + 25 ^1) cos(2cjt + 20o + 6* + 6*1) 

[(f]2 _ ^2)2 ^. Q,2^2]y(^^ 4^2)2 ^ 4^2 jf ' 

<j?)l(pLUj), tan(^i) = (f}2 _ 



where tan(6') = 
^u?)I{2olJ) and 

Aq^Ax-A2 

Ax = WaG^{\) 

A2 = 



WoG2»2 
2 

f72 _ 

1 4" 



•^2)2+^2^21 
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^ WoG2(A)v/(172 _ ^2)2 ^_ ^2q,2 
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Substituting Eqs. ||HJ) and © into Eq. I|6a|) . we obtain 
that the velocity of the kink up to the second order in e 
is given by 



dX 
'dt 



^qWo sin{ujt + 9o + S) h (A) 



i?2/i(A)Va2 +c^2 /^(A) 

-qWi sinjojt + 9q + 5) _ hjX) 
i?2/i(A)Va2 +t^2 /^(A) 



W^2 



.(16) 



By taking the average value of this velocity over one pe- 
riod (T = 27r/w), {dX/dt) = {V) = (l/T) J^{dX/dt)dT, 
we finally obtain 



{V) = 



-e^qWoG{X)n^a 



2i?2/^(A)(a2 ^ ^2)[(J^2 _ ^^2)2 + Q,2^2 



(17) 



Eqs. (|16l) and (|17|l represent the main result of the 
paper. From their analysis the following important con- 
clusions can be made. First, we notice that the non-zero 
average velocity is due to the effective interaction between 
the translational {P{t)) and internal mode (Wi(t)) rep- 
resented by the first term in the second bracket of 
Second, Eq. (|17|) shows that the average velocity does 
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not depend on the initial phase. Indeed, the direction of 
the motion is determined only by G(A), i.e. by the pa- 
rameter A which controls the asymmetry of the potential. 
It is not difficult to check that G(A) > (G(A) < 0) for 
A > (A < 0) and that G(A) vanishes at A = 0. Third, 
from H17|l we also see that for a given frequency there is 
an optimal value of the damping to achieve the maximal 
mean velocity (see Fig. 4 of 14]). This optimal value can 
be easily calculated as a function of lo and A (see the solid 
line of Fig.El), and is given by 



we obtain 



ao,t = y , (18) 

where D = {VL^ - uj'^f . 

Fourth, the activation of the internal mode alone, with- 
out any phonons present in the system, does not allow the 
damping to rectify the motion. By varying a in Eq. H17() 
we can not change the sign of the average velocity. This 
feature confirms the prediction of where the recti- 
fication of the movement for small damping and large 
amplitude of the ac force was related to the excitation of 
phonons for a given choice of parameters. 

Moreover, the existence of a resonant behavior of the 
mean velocity as a function of the frequency is also 
qualitatively confirmed by Eq. (|17|l . As reported in Ref. 
[iJI (see Fig. 6 of this paper), {V) becomes maximum 
when uj approaches the internal mode frequency fJ/ w f2. 
This agreement, however, is only qualitative since the 
resonant peak is very close to the edge of the phonon 
band so that phonons are easily excited in the system. 
In fact, Eq. ((TUl together with Eq. ifTs)) show that the 
kink oscillates with two frequencies to and 2lj. Close to 
the resonance 2lj ~ 2D, one is inside the phonon band so 
that phonons get excited and the CC analysis becomes 
unaccurate. 

Notice that Eq. (|17|l is valid for small e and for times 
(t >> a~^), thus for the zero-damping case nothing can 
be inferred about the mean velocity. It is of interest to 
investigate the zero damping case separately so to check 
the role played by the damping in the phenomenon. For 
a = the momentum equation |(7J simplifies as well as 
the collective coordinate equations (|6a|l - lj6b|l . From Eq. 
O (with a = 0) we have that a solution satisfying the 
initial condition P(0) = is readily obtained as 



Pit) 



9e 



cos(wt + ^o) — cos(6'o)]. 



(19) 



By substituting this equation and Eq. © into Eq. ibai 
we obtain that the kink velocity at order e is 



dX 



'cos{Lut + 9o)-coseo]-^^eWi. (20) 



/i(A) 



The first order correction to the soliton width, Wi, can 
be calculated solving Eq. H1U|) with a = 0. For oj << fl 



WoGiX)cosdoy/Lo^ + nHi,ii^9Q 
Wi{tj = . „ 2^ sm(ilt + tp) 



WoG{X) 



(21) 



where tamp — (fl/ lu) tanOo in order to fulfill the initial 
conditions W^i(O) = and Wi{0) = 0. It is clear that 
the velocity of the kink depends on the initial phase, ^0 
so that, by averaging over the phase, one obtains zero 
transport. A straightforward calculation shows that the 
same result is true also at order e^. This implies that 
the ratchet effect in the zero-damping case does not ex- 
ist. As we will see in the next section, this conclusion is 
confirmed by numerical simulations of Eq. JQ) . 



III. NUMERICAL SIMULATIONS AND 
DISCUSSION 

The CC analysis of the previous section neglects the 
presence of phonons in the system and is based on a par- 
ticular ansatz for the soliton shape in Eq. |(SJ. Moreover, 
the approximated solution H17|) is valid only when the 
perturbations are small enough (e << 1). To check our 
results we compare numerical solutions of Eqs. H()a|) - (|6b|) 
and the approximated solution in Eq. H17|) . with direct 
numerical integrations of Eq. . Numerical simulations 
of (0) were performed by using a 4th order Runge-Kutta 
scheme with steps At = 0.01, Aa; = 0.1 in the finite 
length domain x £ [—50,50], taking into account 200 
time periods. In order to obtain the numerical solutions 
of the CC equations we have integrated Eqs. 16a()- (|6b|) 
with the routine DIVPRK of the IMSL library QS, which 
uses the Runge-Kutta- Verner sixth order method. 

Figures ^ and 13 show the mean velocity dependence 
on the damping coefficient for different values of the am- 
plitude and frequency of the ac force. From Fig. Q we 
see that for low values of e there is a very good agree- 
ment between the CC results and PDE simulations in the 
whole range of a. However, when the amplitude of the ac 
force is increased, the approximated solution (|17|l devi- 
ates from the numerical solution of the CC equations and 
from the PDE simulations, mostly in the optimal damp- 
ing region. In Fig. |21 we show the same dependence 
for an higher value of the driver frequency {lo — 0.4). We 
see that although the approximated solution l(T7|l and the 
numerical solution of the CC equations are quite close, 
there is a discrepancy with the PDE results starting from 
the peak and extending to higher values of a, this indi- 
cating that the CC ansatz is valid only for lj « ujph 
{LUph = 1.2699 for A = -0.5). 

In Fig. 121 we check Eq. H18|l for the optimal value 
c^opti^) of the damping constant for A = —0.5. We can 
see that this value increase up to lu ^ 0.7 and after that 
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FIG. 1: Dependence of the mean kink velocity on the damp- 
ing coefficient for two different values of strength of the ac 
force e = 0.1 (upper curves) and e — 0.05 (lower curves). 
Other parameters are fixed as A = —0.5, oj — 0.1, 6*0 ~ n/2. 
The solid lines represent the approximate solution of the CC 
theory (Eq. 11711 ): whereas circles and pluses refer, respec- 
tively, to PDE numerical simulations of Q and numerical 
solutions of the CC equations I6all - 16 bll . 
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FIG. 2: Mean kink velocity versus the damping coefficient. 
Parameters are A = —0.5, e = 0.1, lu = 0.4, do = 7r/2. The 
solid line represents the approximated solution of the CC the- 
ory (Eq. (llTIl l: whereas circles and pluses are the results of 
the numerical simulations of the PDE Q and the numerical 
solutions of the CC equations l|6afl - l|6bfl . respectively. 



ctopti^) decreases. Again, a good agreement between the 
CC theory and PDE results is found at small oj values. 

By increasing the amplitude e of the driver, reversal 
of current can also occur |Q| . This is shown in Fig. 01 
from which we see that as a is increased the mean veloc- 
ity computed from PDE simulations displays a crossover 
from negative to positive values (i.e. a current rever- 



FIG. 3: Optimal value of the damping as a function of the 
frequency of the ac force for A = —0.5. The solid line repre- 
sents the results obtained from Eq. 1)18^ while the circles refer 
to numerical simulations of Eq. Q . 

sal occurred). This phenomenon is described neither by 
the CC analysis (notice that Eq. ifTTIl predicts a positive 
average velocity for all A < 0), nor by the numerical so- 
lutions of the CC analysis depicted in the figure. This 
agrees with the claim made in Ref. that reversal cur- 
rents do not depend on the internal mode mechanism, but 
rather on the existence of phonons in the system. In or- 
der to show the presence of phonon modes when this phe- 
nomenon appears, we have plotted in Fig.|31the discrete 
Fourier Transform (DFT) of the kink's width W{t) (ob- 
tained from the numerical simulations of the PDE equa- 
tion) for two values of a, i.e. a = 0.2 (before crossover 
occurs) and a = 0.7 (inside the rectified motion region). 
In the former case (upper panel), one of the frequencies 
of the oscillations of W{t) lies inside the phonon band, 
so phonons are clearly excited. In the latter case (lower 
panel) the main frequencies of the spectrum are well be- 
low the lower phonon edge (wp/i/(27r) « 0.2). From this 
we conclude that phonon modes are important for current 
reversals and a theory based on the internal mode alone 
(as the one presented here) cannot describe properly this 
phenomenon. 

In Ref. a resonant behavior of mean velocity as a 
function of ui was also reported (see the peak of {V) at 
w « 0/ in Fig. 3 of this paper). This feature is also con- 
firmed by Eq. ifTTjl (notice that the denominator of this 
expression has a minimum at u; = 17 « f2/) although, in 
this case, the agreement with PDE results is only qual- 
itative. This is shown in Fig. from which we see the 
presence of a resonant structure, with good agreement at 
small frequencies and large deviations from PDE results 
at higher uj > 0.6 values. The agreement at low fre- 
quencies can be understood from the fact that phonons 
in this case are hardly excited and the CC description 
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FIG. 4: Mean kink velocity versus the damping coefficient 
for A = -0.5; e = 0.5, ui = 0.4, 6*0 = 7r/2. The solid line repre- 
sents the approximated solution of the CC theory (Eq. 1)17^ ') 
while circles and pluses denote the results of the numerical 
simulations of the PDE Q and of the numerical solutions of 
the CC equations Iltial - ll6b| l. respectively. The PDE results 
show that for very small values of damping (a < 0.1) pairs of 
kink and antikink appear. 



becomes accurate. On the contrary, when lo gets close 
to the resonant peak, phonons are easily excited and the 
CC analysis becomes unadequate. 

We have also investigated the dependence of the phe- 
nomenon on the asymmetry parameter A as well as the 
importance of the damping term for soliton ratchets. In 
Fig. [7| we show the mean velocity as a function of A for 
fixed system parameters and two different values of a. 
We see that the curve is anti-symmetric around the ori- 
gin meaning that the sign of A determines the direction 
of motion, the maximal effect occurring around |A| — 0.5, 
i.e. the point of maximal asymmetry of the potential. 

Finally, we have investigated the zero damping limit 
of the phenomenon. In FigO we depict the time evolu- 
tion of the center of mass of the kink for a = and for 
different values of the initial condition phase 9q of the ac 
force. We see that X{t) is basically a linear function of 
time, whose slope depends on cos{9q). This implies that, 
in this case, is the initial phase which determines the di- 
rection of the motion and not only the asymmetry of the 
potential, as predicted by the CC analysis. Since in most 
experimental contexts initial phase are usually unknown, 
one must consider the phase as a random variable and 
take average on it. This obviously implies that soliton 
ratchets cannot exist in the zero-damping limit. 
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FIG. 5: DFT of the width of the kink, W{t), obtained from 
the numerical solutions of the PDE for the same parameters 
of Fig. 2] and a = 0.1 (upper panel) and a = 0.7 (lower 
panel, when the current is rectified). The upper panel refers 
to the case in which a current reversal occours. The relevant 
frequencies in the spectrum are: lo\ = 2lo = 0.8; UJ2 = -ii-o = 
1.2; CJ3 = a; = 0.4 and uJi = (7/2)aj = 1.4 (this last being 
inside the phonon band). The lower panel correspond to the 
case in which the current is rectified. The main frequencies 
in this case are located at cji = w = 0.4 and a;2 = 2tj = 0.8 
away from the phonon's band. The frequency of the internal 
mode and the lower phonon's frequency are Q,i — 1.056 and 
i^ph ~ 1.269, respectively. 



IV. CONCLUSIONS 

In this paper we have studied the ratchet dynamics of 
the kink solution of the ADSGE by using a collective co- 
ordinate approach with two collective variables, the cen- 
ter of mass and the width of the kink. For these variables 
we have obtained a system of ODEs from which we de- 
rived an approximated expression for the mean velocity of 
the kink as a function of the system parameters. We have 
confirmed that the CC approach is valid when phonons 
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FIG. 6: Mean velocity of the kink as a function of u. 
The pluses superimposing the solid line show a good agree- 
ment between the approximated (solid line) and the numeri- 
cal (pluses) solutions of the CC theory. Results obtained from 
the integration of the Eq. Q (circles) coincide with the CC 
theory only for smaller values of uj. Parameters are fixed as 
A = -0.5, e = 0.1, 6*0 = 7r/2 and a = 0.5. 
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FIG. 7: Mean velocity of the kink center of mass versus A 
for e = 0.1, Lu = 0.1, 6o — 7r/2 and for two values of the 
damping coefficient: a = 0.2 and a = 0.5 (for fixed A, the 
mean velocity of the a = 0.5 case is smaller). 

are not excited in the system, i.e. for small values of e 
and for w << ujph- We have shown that for a proper de- 
scription of soliton ratchets it is not enough to consider 
the kink as a point particle moving in a ratchet potential 
0,113, but it is crucial to include also the internal oscil- 
lations of the kink profile. In particular, we have shown 
that the net motion of the kink becomes possible when its 
internal and translational modes become effectively cou- 
pled (the effective coupling being possible only in pres- 
ence of damping). We also showed that the asymmetry 
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FIG. 8: Time evolution of the center of the kink for different 
values of initial phase, So, in the non damped case. Parame- 
ters are fixed as: A = —0.5; e = 0.1 and lo — 0.4. The solid 
and the dashed lines represent the numerical simulations of 
Eq. Q and the numerical solutions of Eq.s I6all - ll6b^ . respec- 
tively. 

of the potential determines the direction of the motion 
and that in the zero-damping case the ratchet effect van- 
ishes (i.e. it depends on initial conditions). The resonant 
behavior of the velocity as a function of frequency and 
damping was also investigated. We found that the mean 
velocity approach a maximum value when the frequency 
of the ac force goes to the internal mode frequency or 
when the damping coefficient approach its optimal value. 
Finally, we have shown that the occurrence of current re- 
versal is related to the presence of phonons in the system 
rather than the coupling between the translational and 
the internal mode. 

In conclusion, the results of our analysis confirm the 
internal mode mechanism for soliton ratchet proposed in 
Ref. Ql and provide an approximate analytical descrip- 
tion of the phenomenon. 
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